Self-similar dynamics of bacterial chemotaxis 
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We investigate the pattern formation of colony generated by chemotactic bacteria through a 
continuum model. In a simplified case, the dynamics of system is governed by a density-dependent 
convection-reaction-diffusion equation, u t = (u m ) xx — 2n{u m ') x + u — u m . This equation admits 
the analytical solutions that show the self-similarity of the bacterial colony's morphogenesis. In 
addition, we found that the colony evolves long time as the sharp traveling wave. The roles of 
chemotaxis on the regulation of pattern formation in these results are also discussed. 

PACS numbers: 87.18.Hf, 05.45.-a, 87.23.Cc, 82.40.Ck 



The spatiotemporal pattern formation in bacterial 
colony is an example of elegant self-organization in na- 
ture. This phenomenon attracts much attentions of sci- 
entist in interdisciplinary fields for decades [10. The 
pattern formations are not only resulted from reproduc- 
tion and movement in bacterial populations but also from 
the cooperation of the individual to cope with the hos- 
tile environmental conditions. This collective behavior 
reflects the bacterial communication capabilities and so- 
cial intelligence which might promise the newparadigm 
for biocomplexity and natural computations 

(EH- The 

underlying mechanism of bacterial pattern formation is 
a key to understand living organisms, which is therefore 
important to biotechnology and medical science. 

It has demonstrated that the colony of various bac- 
terial species exhibits branching pattern P, [H-Q which 
is analogous to the fractal pattern formation in the 
diffusion-limited aggregation (DLA) process Q. In this 
manner, the bacterial colony's patterns are typically self- 
similar, where they are the same at every scale (scale 
invariant). This macroscopic pattern is arisen from the 
microscopic behavior of each individual member of bac- 
teria in colony. The bacteria move in fluid medium by 
swimming as random walk motion, in which the bacte- 
ria propel themselves in nearly straight run separated by 
brief tumble. Unlike nonliving particles, the bacteria re- 
spond to a chemical attractant such as nutrient by swim- 
ming along its gradient, known as the chemotaxis. The 
chemotaxis has an essential role on the regulation of bac- 
terial colony pattern formation as demonstrated by Refs. 
[§- 11 1 . The bacteria detect the spatial gradients by com- 
paring a temporal difference between the amount of at- 
tractant molecules that bind to the membrane receptors 
along their path. The bacteria implement the chemo- 
taxis by delaying the tumbling frequency as cells swim 
up the gradient of the attractant or down the gradient of 
repellent. Moreover, the bacteria perform the intricate 
cellular communications capacities in forming the colony 
to respond the environmental conditions: These includes 
long range chemical signaling, such as quorum sensing 



and chemotactic signaling. At first glance, the formation 
of bacterial colonies appears quite complicated; but the 
reaction-diffusion model has successfully described the 
continuum limit of this dynamics 0, 0, H 0| ■ 

In this Letter, we consider a nonlinear reaction- 
diffusion model that has been proposed for studying 
the pattern formation in bacterial colonies exemplified 
by Paenibacillus dendritiformis grown on Petri dish 
[E 0, E2, EH ■ This bacteria species is motile on the dry 
surface by cooperatively producing a layer of lubrication 
fluid in which they swim and its colony exhibits branch- 
ing pattern. Indeed, the numerical simulations on this 
model has been carried out from past studies 0, EH, EH ■ 
It seems likely that the numerical results reproduce the 
branching pattern in the bacterial colonies, which agrees 
well in comparison with experimental data. However, the 
sole numerical results are inadequate to explain the scal- 
ing laws of bacterial colony growth. Thus, an analytical 
result is still required to archive this task. The aim of this 
work is to analyze the simplified form of those proposed 
nonlinear reaction-diffusion models for pattern formation 
in bacterial colony, which takes chemotaxis into account 
[E S EH E3| • Our analytical results could be plausible 
for interpreting the results in both experiments and sim- 
ulations. 

We now explain the bacterial chemotaxis model in our 
consideration. Actually, the bacterial colony evolves in 
two dimensions however; each tip grows in one dimension, 
except for occasional branching. Thus, it allows us to in- 
vestigate this problem in one-dimensional space which its 
results could be equivalent to one obtained from the two- 
dimensional system Q. As proposed in Ref. 14 1, the 
dynamics of bacterial populations is governed by a gen- 
eralized convection-reaction-diffusion (CRD) equation 
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where u(x,t) and s(x,t) are, respectively, the bacterial 
density and the attractant density in spatial coordinate 
x and time t. D(u), R(u) and i9(s) are, respectively, the 
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diffusion coefficient, the reaction term and the drift ve- 
locity due to the chemotaxis. This equation is similar 
to the generalized Keller- Segel equation tl a. E q. is a 
simplified form of the full model in Refs. [7|, [l2j, ll3| under 
following assumptions, i) The nutrient density is propor- 
tional to the bacterial density and is absorbed some way 
in reaction term, ii) The production of lubricant fluid 
is proportional to the bacterial density and is absorbed 
into the medium. Thus, the lubrication fluid can be re- 
moved from the dynamics and it is represented through 
the diffusion coefficient 12|, |f4[ . iii) The chemotactic sig- 
nal can be also a field produced directly or indirectly by 
the bacterial cells. 

In the usual case, the diffusion coefficient is constant 
which represents the Brownian motion or random walk. 
Since the random walk seems to be unrealistic move- 
ment for the active individual such as members of bi- 
ological population, the crowd-avoidance movement has 
been introduced [1 d - fecf . In this case, the diffusion 
coefficient is represented by a density-dependent form 
D(u) — M(u/a) p , where M > is diffusion constant, 
<j = limt^oo u(x, t) is equilibrium density and p > 0. For 
the reaction term, it can be represented by the general- 
ized logistic law R(u) = au[l — (u/a) p ], where a > 
is rate constant [20j. The chemotatic drift velocity can 
be expressed as 'd(s) = C( s )x( s ) s x, where x( s ) s x acts as 
the gradient sensed by the bacterium (with -y(s) having 
the units of 1 over the chemical concentration) [12| • C( s ) 
is the bacterial response to the sensed gradient and has 
the same units as a diffusion coefficient [f2j]. Therefore, 
we assume £(s) = 7-D(it) = r yM(u/a) p where 7 is a con- 
stant, positive for attractive chemotaxis and negative for 
repulsive chemotaxis [f2[ • Here, we are interest in a spe- 
cial case where the gradient sensed by the bacterium is 
nearly uniform and x( s ) s x is treated as a constant [Til ]. 

By substituting D(u), R(u) and i9(s) into Eq. (JT]) with 
the transformations t* = at, x* — (ma/M)2x, u* = u/a 
and k = (l/2)7(mM/a)5^(s)si, we obtain the dimen- 
sionless equation 



u t = {u m ) xx - 2n(u m ) x + u-u n 



(2) 



where m = p + 1 > 1 and the asterisk is dropped. 
The traveling wave ansatz to Eq. @ has been analyzed 
14i |2I| . It found that there exists traveling wave solution 
for wave speeds equal to or greater than the minimum 
value [1J, [21| . Although the solution of Eq. ([2]) has well 
understood as the traveling wave, the exact or explicit 
solution in space-time coordinate has been unknown. In 
addition, it has been pointed that the traveling wave so- 
lution is only intermediate asymptotic behavior of system 
[22I]. Therefore, it is still unclear that how the systems 
evolve to traveling wave from an initial state. In previ- 
ous study by Ref. [23], Eq. ((2]) can be mapped to an 
inhomogeneous diffusion equation by a transformation 
technique, where the reaction term is eliminated. How- 
ever, the exact solution has not been obtained in the 



closed- form. In our previous work 2J], Eq. ([2]) without 
chemotaxis can be differently mapped to an anomalous 
diffusion equation [25I - I27I ] which the exact self-similar so- 
lution can be obtained for any m > 2. We then extend 
the similar technique to analyze Eq. ©. 

We perform the analysis on Eq. © with similar 
method to Refs. [HIH. Note that Eq. ^ can be writ- 
ten in equivalent form as 



d__,\ _J_fd_ 
dt J U ^ \dy 



(3) 
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where y — x/uj and uj = h 
Next, Eq. (|3]) can be evaluated to 

To eliminate the ex- 
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ponential terms in Eq. ([3]), we define the solution 

as 



u(y,t 

By substituting Eq. 
e -(™-i)*$ t 



e*e 



(4) 



»*(y,t). 
into Eq. ©, 



(4) 

we obtain 



[e {u +1)y (t> m ) y }y By intro- 
ducing the variable transforms dr = (m — 1 )e( m ~ 1 ^ t dt 
and d(f> = ^+^L e {rn+^ 2 )y/rn dy that giveg 



r{t) = e ( m - 1 )* 



f 



this equation can be reduced further to 



where k 



1 ( m+ui 2 V 
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Eq. © has the form of anomalous diffusion equation 



as studied in Refs. |25l - l27| . To obtain the solution of 
Eq. ©, we assume the solution to be the scaling func- 
tion of the type 



-F 



F{9) 



(7) 



T(r) ! 
By substituting 



-^(F + 9F e ) 



T(t)~ \T(t) 

where 6»(</>,r) = <j)/T(r) @-[27| 
Eq. © into Eq. ([5]), we obtain 
T Jl l+2 \9 l (F m )o\ which can be separated into two equa- 
tions: T m - l T T = A and -\(6F) 9 = k [0 l (F m ) e ] e where 
A is constant. Performing calculations similar to Refs. 
2J-|27[ and choosing A = (m — I + l)^ 1 , we obtain 



T(t) = (t + a ) ™('"-d 
F{6) 



(8) 



where a > and b are constant. By substituting Eq. ([5]) 
back into the main solution Eq. (j4]), after doing some 
algebra, we obtain the explicit solution to Eq. ([2]) 



u(x, t) 



(e pt - 1 + a) * 



pu)x 



(e pt - 1 + a) u 



(9) 
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We now consider the initial condition for solution (|9]) 
uq(x) = u(x,0) that has (semi-jcompact support: 
lim^^-oo uq(x) = p and Uq(x) = for x Xq, where 
p is initial density amplitude and Xq is initial front po- 
sition. Such conditions, we obtain a = p~ p and b — 
_ p -pu 2 /(p+i) e -pux /(p+i) _ By substituting the constant 
a and b into Eq. ([9]), we obtain the exact solution to 
Eq. © 



u(x, t) 



X il 



pe 



[pp( e p* - i) + i]p 

e pu){x-x ) 

[pp( e pt - 1) + l]" 2 



(10) 



Since the solution Eq. ([10)) has two forms depending 
on the value of w, we define u+(x, i) and u-{x, t) are the 
solutions corresponding to ui+ = k + y/ K 2 + 1 and w_ = 
k — \/k 2 + 1, respectively. We now construct another 
solution of Eq. © by linear combination of these two 
solutions: 



uu(x, t) = it+(x, t) + U-(x, t). 



(11) 



By using binomial expansion, we obtain an approxima- 
tion 



0(H), (12) 



where H = Y^=i (™)u k +u m ~ k . By substituting Eq. CP 
as solution into Eq. ([2]) and using Eq. (fT2|) . we obtain 
w t = K ra ) :c:c -2i/(u; m ) :c + W - W m + 0( J ff :c:!; -2 J /H ;c -//). 
Therefore, io(:r, t) in Eq. (|TT|) can be a solution of Eq. ([2]) 
if the correction 0(H XX — 2vH x — H) can be ignored. 
Using Eq. (fT2"j) . we obtain the solution 



w(x,i) « 2~p [u^.(x,t) + i£(:r,t)] : 



(13) 



where (2) 1 / p is normalized factor. We note that in case 
of no chemotaxis k — 0, thus a; = ±1, these results have 



been discussed in our previous work [24 [ . 

The evolution in space and time of bacterial density 
profile U-\-{x,t) and U—(x,t), as in Eq. (|10[) . are illus- 
trated in Fig. (P). It is seen that the density profiles 
start from the initial state uq(x) then grow and expand 
to the unoccupied region. At sufficient large time scale, 
the density profiles reach the saturated point with the 
maximum value 1. After that, they seem to propagate 
with unchanged shape; u+(x,t) is propagating to the 
right whereas it_(x, t) is propagating to the left. Roles of 
chemotaxis on the regulation of pattern formation in the 
system is reflected by parameter k. Since k < V k 2 + 1, 
cj + is always positive whereas u;_ is always negative that 
leads the tails of u+ to decay as x — > oo and the tails 
of u- decay as x — > — oo. In the former case, the front- 
interface is sharper because |w+| > |w— |. Due to the in- 
fluence of chemotaxis, the distribution of density profile 
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FIG. 1. The spatiotemporal evolution of the density profile 
u(x,t) (Eq. (JTOj)) in the case of p = 2 with initial conditions 
p — 0.2 and a;o = 1. The solid lines represent u+(x, t) and the 
dashed lines represent U-(x,t). 




FIG. 2. The spatiotemporal evolution of the pulse-like density 
profile w(x,t) (Eq. (|13[0 in the case of p — 2 with initial 
conditions p = 0.2 and xo = 1. 



is biased toward to the right thus the front of w+ is mov- 
ing faster than of u_. The spatiotemporal evolution of 
the combined density profiles w(x,t) is also illustrated in 
Fig. ([2]). The density w(x,t) form the pulse-like profile 
that grow and expand with asymmetric shape. It be- 
haves like u + (x,t) for x 3> xq and behaves like u—(x,t) 
for x <C — xq. Due to the bias force from chemotaxis, the 
peak of w(x, t) is moving toward to the right. 

From Eq. (fTU|) , we calculate the front position r(t), 
that the density falls to zero u(r, t) = as r(t) — xq + 
uj ln ^ p ^ e . The plot of relative front position r(i) — 



Xq is shown in Fig. (J3]). The relative front position of 
U—(x, t) is slow varying when compared with of u+{x,t). 
At sufficient large time, that e pt ~3> 1 and (Pe vt ^> 1 
thus t' w — In p, the relative front position seems to vary 
linearly in time r(t) — xq ~ ut. It implies the constant 
front propagating speed. Consequently, we calculate the 
front speed as v(t) = f t r(t) = pP "Jtl"i )+1 - At large 
time scale t 3> t', the spreading speed trends to be the 
constant c = lim^oo v(t) = lu(k). At this point, it is 
clearly seen that the spreading speed is biased by the 
chemotaxis through the parameter k. 

At the large time scale, t ^> t', the bacterial density 
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is nearly uniform. The solutions show explicitly the self- 
similar structure of the bacterial colony growth. Con- 
sequently, the scaling law of bacterial colony growth is 
reveled. As usual, the solutions evolve long time scale as 
the traveling wave. 

K. Khompurngson acknowledges the Centre of Excel- 
lence in Mathematics (Thailand) for partial financial sup- 
port. 



FIG. 3. The relative front position r(t) — xo corresponding to 
density profile in Fig. ([T]). The solid lines represent the rela- 
tive front position of u+(x,t) and the dashed lines represent 
the relative front position of U—(x,t). 



profile Eq. (|I0[) emerges the traveling wave form 



u(x — cut) 



S(i-ui-io) 



ep+ 



(14) 



where uj is front speed. When converting to the physical 
dimension, the front speed obtained here is comparable 
to the minimum value c m in — Wi/ 1^ for the sharp trav- 



eling wave front [LJ, l2l| . Similarly, at the large time scale 
t 3> t', Eq. (fl~3"|) develops to the expanding pulse-like wave 



w(x — U)±t) 



_(x — ijJ+t) + vP_{x 



-*)]' 



(15) 

with the expanding speed oj±. 

We point out that, at large time scale t 3> t\ the so- 
lution Eq. forms a scaling law $(0, r) « -^F(^), 

I 2 

where B = , . In biological context, the solution 

" m{m — 1) ° 1 

captures the self-similar features of the spatiotemporal 
pattern formation in bacterial colonies in the terms of 
transformed quantities: u — > e" x / m e*$, r — > e (™- 1 )* ; 
and (j> e (™+" 2 )(*-*o)/m^ This se if. s i m ii ar solution 

converges to the traveling wave solution which its behav- 
ior can be classified as the intermediated asymptotics of 
the second type, as described in Ref. p2j ]. 

In summary, we investigated the spatiotemporal pat- 
tern formation of bacterial chemotaxis described by 
a density-dependent convection-reaction-diffusion equa- 
tion. The analytical solutions are obtained for the sim- 
plified case, where the gradient sensed by the bacterium 
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